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Within the frameworks of spin-polaron concept and the spherically symmetric state 
for the antiferromagnetic spin background, the one-particle motion is studied for two- 
dimensional Kondo-lattice. The elemetary excitations are represented as a Bloch super- 
position of four one-site electron states: two local states- a bare electron state and a 
local spin-polaron of small radius, and two states of delocalized polarons which corre- 
spond to the coupling of local states to the antiferromagnetic spin wave with momentum 
Q — (n, 7r), so called Q-polarons. As a remarkable result we show that the lowest band of 
elementary excitations is essentially determined by Q-polaron states in strongly coupled 
regime. The account of Q-polarons shifts the band bottom from (n,n) to (0,0). The 
spectral weight of a bare particle in the lowest band states can greatly differ from 1. This 
may lead to a large Fermi surface for relatively small particle concentration. 
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In order to understand the nature of high temperature superconductors it is important to describe 
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properly the motion of a hole in the CuOi plane . This motion takes place on the antiferromagnetic 
spin background of copper spins and it must be treated as a correlated motion of a hole coupled to 
spin excitations (a spin polaron) . Usually the spin polaron is studied within the frameworks of the t- J 
model |l| and the three-band Hubbard model These models seem to be rather complex because 

the effective spin-fermion coupling appears there only through the fermion hopping. Moreover, in the 
case of t-J model it is difficult to take into account the local constraint which represents the strong on- 
• site fermion repulsion. We shall study the formation of the spin polaron within the frameworks of the 
I more simple Kondo-lattice model. Note, that a generalized variant of the model leads to description of 
I/"") i low- lying states in the three-band Hubbard model . 

Treating this model Schrieffer S replaced the spin subsystem by a classic field with a doubled lattice 
period and described the elementary excitations as the superposition of electronic states with the momenta 
p and p + Q, where Q = (tt, 7t) is the antiferromagnetic vector. Such a description is deficient for the two- 
dimensional S = j antiferromagnet because of strong spin fluctuations. Moreover, it becomes impossible 

S if the spin subsystem is found in the homogeneous spherically symmetric state with long range order 
but with zero value of the average site spin < S a >= (the homogeneous Neel state) [p|-[ll[. The 
■ homogeneous Neel state seems more adequate, then the usual two-sublattice Neel state, for treating the 
£^ ' ground state of the doped CuO-i plane. 

The distinctive feature of the present investigation consists in considering of the one-electron motion 
^ ' on the background of the homogeneous Neel state of the Kondo lattice. This motion is described by a 
^ , spin polaron with the spectrum periodicity relative to the full Brillouin zone. Note that the conventional 
two-sublattice spin approach leads to periodicity relative to the magnetic Brillouin zone [^||] . 

Another distinctive feature of our investigation consists in treating the spin-polaron as a complex 
quasiparticle accounted for a coherent superposition of a bare electron, of a local polaron and of two 
types of antiferromagnetic delocalized polarons. A local polaron can be attributed to the local singlet 
(an analogous of Zang-Rice polaron in the three-band Hubbard model) and such a polaron represents 
the lowest band in the limit of the strong Kondo interaction J. As to the antiferromagnetic delocalized 
polaron it is a bound state of an electron (or of a local polaron) with a spin wave with the momentum 
Q. As we know such bound states were not investigated previously. If the spin subsystem is found in the 
state with the antiferromagnetic long range order then the amplitude Sq of the spin wave with q = Q 
(the Q-wave) has the macroscopic large value and has the properties analogous to the amplitude of a 
Bose particle with the zero momentum in the supcrfluid Bose-gas. As a result for many problems this 
amplitude can be treated as a c— number. Then the coupling of the Q-wave to local electron states does 
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not represent new states but leads, as mentioned above, to mixing of the states with the momenta p and 
p + Q. But this treatment fails in the case of the homogeneous Neel state. In this background the average 
value < Sq >= and only < SqSq > can be treated as a macroscopic value. Then the coupling of a 
particle local state to Sq corresponds to a new delocalized state. It will be shown that it is important to 
take into account such a quantum nature of the spin Q-wave because the transitions between the local 
polaron states and the delocalized polaron states essentially determine the spin polaron bands. 
The Kondo-lattice Hamiltonian has the form 

H tot = H + H 1 +H 2 , (1) 
H o = Yl tgat+gUr = e P a+a p , (2) 
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Hi = J^^a+ S r a r , (3) 

r 

H 2 = llJ2 S ?+9 S ?- W 
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Here the sums run over the sites r of a square lattice and over the nearest neighbors with the lattice 
spacing | g |= 1. For short we miss the spin index for creation a^ a and annihilation a ra operators of the 
Fermi particles (we shall name them as electrons) and in the Hamiltonian of Kondo-interaction H\ we 
take the notation S r = S"a a . 

Let us represent the first two equations of the infinite chain of equations for the retarded Green's 
functions, which describe the motion of one particle in the antiferromagnetic spin background. In the 
present paper we restrict ourself mainly to the analysis of the spectrum of the quasiparticles. That is why 
for simplicity of notations we shall simply write operators a instead of the Green's functions << a \ b + >> 
and shall miss the nonhomogeneous terms < [a, + >. The equations have the form 
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Jb p , b p = AT* b r e- ipr , b r = S r a r , (5) 



uiS r+ na r — [S r +Ra r , H + Hi + H 2 ] = 

= tgS r +Ra r+g + J S r+ ftS r a r + He a fjj a a^r+R+g^r+R ar - (6) 

g g 

In Ref. H the mean-field approach was taken for the Neel lattice and the spins were treated as the 
classical vectors: 

5 r Q = S az S Q e lQr , So = const, Q = (±tt, ±tt). (7) 

In this approximation the Kondo-interaction Hamiltonian (Q) takes the form of the potential energy 
with the doublet period and the equation (^) has a closed form: 

(u) - e p )a p — aJS a p -Q. (8) 

Here p is in the Brillouin zone with the periods (2tt, 0), (0, 2ir), but the spectrum has the periodicity of 
the magnetic Brillouin zone: 



E ap = ±yje* + A 2 , A = S J. (9) 
In (|^) and below we suppose the spectrum model with the "nesting": 

e p = ~ £ p+Q = —2t(cosp x + cospy) (10) 

In the case of the S — \ spin system the quantum fluctuations are important and the equation (^|) 
must be used in order to find the last term in Eq.(||). In one's turn the equation (^) hasn't a closed form. 
In order to close the chain of equations we use the standard Mori-Zwanzig projection technique fl2]| . In 
our case this means that wc must approximate the last two terms in the right-hand side of (^|) by their 
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projections on the restricted space of bases operators. The choice of the restricted set of basis operators 
must be dictated by the physics of the problem under discussion. 

For a paramagnetic spin subsystem state the simplest set of basis operators may be taken as two 
operators that appear in the first equation (|^), that are the operator of a "bare" electron a r and the 
annihilation operator of an one-site spin polaron 6 r .The local polaron excitations in the Kondo lattice 
were previously discussed in Ref. 

Now let the spin subsystem to be found in the homogeneous Neel state with the tensor parameter of 
the long range order: 

< S? >= 0, (11) 



Cf =< S?S? +R >= M a0 e lQR , | R |» 1. (12) 

Then it will be consecutively to introduce the operators which take into account the correlation of the 
electron and the local polaron with the antiferromagnetic Q-w&ve 

c r = Q r a r , c p — Qoa p+ Q, (13) 
Qr = N- 1 e tQR S r+R = e^ R Q . (14) 

R 

d r = Q r S r a r , dp = Qob p+ Q. (15) 

Note, that the operators c r and d r allow for distant correlations between the spin and the electron 
(| R \» 1). So, in a sense our approach is alternative to the widely used t-J model investigations based 
on the decoupling of an on-site Fermi operator into a spinless fermion and an antiferromagnetic magnon 
operator ||^1§. 

It is substantial that the adopted set of the basis operators is closed relative to the Kondo-interaction 
Hamiltonian: 

[a r ,i?i] = Jb r , 
3 

[b r ,Hi] = J(-a r - b r ) , 

[c r ,H 1 ]=Jd r , (16) 
3 

[dr,-ffi] = J{^c r - d r ) . 

Taking into account Eq.(^J) this gives, in particular, 

(cj - € p+ q)c p = Jd p . (17) 

In order to get a closed form for the Green's functions equations we shall project the corresponding 
commutators on the following orthonormal set of basis operators: 

B\ r = a r , 

B 2r = —=b' r , b' r — b r — c r , 

V/2 

h = < [b' r , b'+] + >= | - M, M = M aa , 
B 3r = -^=c r , (18) 

V J3 

f 3 = < [c r ,c+]+ >=M, 

1 M 
B ir = —r=d' r , d' r — d r - Ma r + —(b r - c r ) , 

V/4 J 2 

M 

.f4 = <[d' r ,d'+] + >=M(f 2 -—). 

h 
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The coefficients of the projection relations 



[Bi,H] =J2 a v B i 



(19) 



are determined as 



an =< [[B U H],B+] > 



(20) 



Taking into account the equations (Q) , (^) and ( p0| ) , we obtain in momentum representation a simple 
system of four equations 



(21) 



The matrix elements of the spectral matrix aij are equal 



a(l,l) 



a(l,2) = JV/2, a{l,3) = JVM, o(l,4) = 0, 
o(2, 2) = i- [(C + Af)e p - J(/ 2 - M) - 4JC] , 

J2 



o(2,3) 



M 1 
i(2, 4) = [-{C + M)e p J/ 4 - 27M(1 + — )] -_, 

72 J2 VJ2/4 




o(3,3) 



a(3,4) = J 




o ,n ^r/^ m m 3J IM { 8M M, 4M n 

a(4,4) = _[(C + M)-^ - (C+ — )] £p - — + — [— -4C(1 + ^) - — ]. 
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Here C = C^" 1 < 0. To obtain these expressions we used the following approximation of Takahashi 
.7J for four different-site spin correlation functions 



< m C3 cfe ql > = 



p%3 sykl 
L/ n-r2 U r3-r4 



T*l — T3 — ^4 



(23) 



Hence, in the adopted approximation the electron motion in the antiferromagnetic background is 
described by the quasiparticle that is a coherent superposition of four Fermi fields - the field of the bare 
electron a p , the field of the delocalized polaron c p and two fields of the localized polaron b p and d p , which 
are hybridized mainly due to the Kondo-interaction. The system of equations (|l]) can be treated as 
the system of the Schroedinger equations where the matrix elements aij reproduce the amplitude of the 
transitions from the state j to the state i. The presence of the nondiagonal matrix elements underlines 
the quantum nature of the spin S and the Q-waves. The eigenfunctions and the eigenvalues of Eqs. 
( pl| ) describe the elementary excitations which generate four bands. The form of the bands depends on 
the state of the magnetic subsystem (the quantities M and C) and the relations between the energetic 
parameters t,J,I. Below we take the following typical values for the quantities M and C: M = 0.1 and 
C = -0.335 

In the adopted model all the spectra depend on momentum only through e p . Let us mention that 
calculating the matrix elements (|20| ) we took the approximation of low density of Fermi particles, i.e., 
p >-> 0. 

It is essential that the resulting four bands are only approximately symmetric relative to the magnetic 
Brillouin zone boundary. In the momentum space the distance between two quasiparticle states with the 
equal energy can differ from the antiferromagnetic vector Q. This means that dictated by these bands 
the "shadow band" effect can be other than the usual one displaced by Q vector in the two-sublattice 
Neel antiferromagnet [^9[ . 

Four waves of the quasiparticle spectrum E^\p),i = 1 — 4, given by the eigenvalues of the system 
of equations (|2l]), have a clear physical meaning in the cases when one of energetic parameters t, J, I is 
much greater then the others. The bands for these cases are represented in Figs. l-3(a). In Figs. T3(b) 
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we also represent the values of the electron Green's function residues Z p 



(i) 



which correspond to the poles 
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4 z (i) 

« a pa | o+ »=« Si, P | »= ]T u _E(f){pY (24) 

In the low density limit the residues Z p 1 ^ are determined by the solution of the equations of motion 
for Green's functions << Bi_ p \ B^ p >> when the nonhomogeneous terms are equal to Wi.\ =< 

[Bi tP , Bi p ]+ >= Si. i- The Z p ^ values represent one-particle spectral function A(p,u>) = J2i Z p ^8\u — 

E^(p)]. Each of these values Zp characterizes the contribution (weight) of a bare particle state a p to 

the quasiparticle state with the energy E^(p) . For this reason ^ Z p ^ = 1. 

First, if the energy of the antiferromagnetic exchange is the largest energetic parameter (I » J,t) 
then the formation of the local polaron b r is energetically disadvantageous because this formation leads to 
the disintegrating of antiferromagnetic correlations between spins S r and S r + g . In Fig.[l] the quasiparticle 
spectra E^(p) are given along the symmetry line p ~ p x = p y for I = 5, J = t = 1. The local polaron 
bands (they represent the motions of b r and d r ) lie essentially higher by energy then two other bands which 
correspond to the hybridized spectra (|j) with the gap A = J\f~M. Two higher bands are characterized by 
a small weight of spectral function A(p, lo) (see Fig. 0(b)). As to two lowest bands, the spectral function 
weight is close to 1 when the band represents mainly the motion of a bare particle and is close to zero 
for the Q-polaron band. In the region of strong hybridization, p close to 7r/2, the spectral weight is of 
the same order for both bands. It may be seen, that these two lowest bands have the symmetry close to 
the symmetry given by the magnetic Brillouin zone, but the corresponding quasiparticles have the full 
Brillouin zone symmetry . 

Now let us discuss the second case, represented in Fig.|| by the spectrum for J = J = O.li, when 
the spin interactions are small relative to the kinetic energy of the electron, t >> J, I. In this case all 
quasiparticles are approximately independent: the free motion of the bare electron B\ p — a p is described 
by the band E\ = e p and the antiferromagnetic polaron B^ p moves in the band E% = £ p +q = — e p . 
As to the local polaron Bi p it has it's own dispersion E^ = a(2,2), The absolute value of the 

multiplier (— | C\ \ +M) / fa of the term with e p in the expression for o(2, 2) characterizes the essential 
band narrowing, i.e. the increase of the local polaron mass relative to the mass of the bare electron. 
Finally, the B^ p polaron band E^ = a(4, 4) also demonstrates a strong, but different renormalization of 
the mass. The difference in the renormalization is one of the reasons why the resulting hybridized spectra 
are asymmetrical relative to the magnetic Brillouin zone boundary. The crossing of the nonhybridized 
spectra takes place close to this boundary where e p = 0. Out of this region one of the bands explicitly 
represents the motion of the bare particle and has the spectral weight Z p close to 1. 

Fig. 3 describes the case J >> t, I. The lowest bands (with the center at E^ = — |j) correspond to 
the motion of the local polaron in a singlet state =b—\a and to the motion of such a polaron coupled 
to the Q-wave. The upper bands (centered at E^' = \J) describe the motion of the local polaron in a 
triplet state b^ = b+ |a and also its coupling to the Q-wave. The "singlet" and "triplet" terminology 

suppose, that acting on the spin subsystem state, the operators br + , fe^^create accordingly a singlet or 
triplet spin-electron pair in the site r. 

The qualitative and important difference of this case from the previous ones is the following. For two 
first cases the bottom of the lowest band is determined mainly either by the motion of the bare particle 
or by the motion of Q-polaron c p , the residues being close to 1 and zero correspondingly. In the present 
case the bottom of the spectrum of the elementary excitations turns out to be substantially determined 
by the b p and d p -states. At the same time, as it may be seen from Fig.||(b) (the solid line), the possible 
filling of each p state by a particle with a fixed spin is not too small for the lowest band and is close to 
0.1. Let us mind that such a strong deviation of filling from 1 must lead to relatively great Fermi surface 
even at small filling. 

The present case demonstrates the importance of taking into account the Q-polaron d p . In order to 
illustrate this circumstance, in Fig. [I| we show the spectrum of the elementary excitations calculated 
for the same values of energetic parameters as in Fig. 3, but using the basis of the first three operators 
a p , b p , c p . As it may be seen by comparison of Figs. 3 and 4, this leads to the shift of the lowest band 
bottom from the point p = (0, 0) to the point p = (ir, ir). At the same time the value of the residue Z p 
for the band bottom is close to 0.1 as in Fig. 3. 

Finally, if three energetic parameters are of the same order of magnitude then the elementary excitation 
in any of four bands is a coherent superposition of the states a,b,c,d and these states have relatively the 
same weight. It is notable that the presence of the long range order doesn't lead to the periodicity of the 
spectrum relative to the magnetic Brillouin zone. 
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Let us mind that the described above properties of the spin polaron are mainly preserved if we suppose 
that the spin subsystem has no long range order but the spin correlation length is large. This must lead 
to the replacement of (|l^) with the following relation which is true for 2D antiferromagnet at T ^ 0: 

Cf =< S?S p r+R >= M^eWe-V, L » 1. 
In this case the delocalized polaron is described by the relation 

c r = L~ 3 V" e lQR e~^ S r+R a r . 
R 

In conclusion we want to make several notes. First, above we investigated the elementary excitations 
of the system - one electron plus the AFM spin subsystem. It may be seen that all the results for the 
system with one hole may be obtained from the discussed above case by substituting (— t g ) instead of t g 
in the kinetic Hamiltonian. 

Secondly, the ignoring of the Q-polarons, i.e., the states c and d, leads to the disappearance of bands, 
and can change essentially the remaining bands. We think that it is important to take into account 
Q-polarons in such models as t— J and three-band Hubbard model, where the particle motion is strongly 
coupled to the spin subsystem. 

Finally, let us mention once more that our calculations are performed in the limit of low filling. Nev- 
ertheless the calculated upper bands have physical meaning if we shall study the electron transitions. In 
general case it may be seen that the one-particle spectral density strongly depends on filling. 
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FIG. 1. (a)The quasipartical spectra E^\p) iorI = 5^>J = t= l along the symmetry line p = p x = p y ; (b) 
the values of the electron Green's function << a pa | a pa >> residues Z p ^ which correspond to the poles E^\p). 
Different line types represent four bands and corresponding residues. 

FIG. 2. (a)The quasipartical spectra E^'(p) for t = 1 3> / = J = 0.1 along the symmetry line p = p x = p y ; (b) 
the values of the electron Green's function << a pa \ a pa >> residues Z p ^ which correspond to the poles E^\p), 
Zp ^ are shown for the interval close to p — tt/2. Different line types represent four bands and corresponding 
residues. 

FIG. 3. (a)The quasipartical spectra E^(p) for J = 5 S> t, I,t = 1,1 = 0.1 along the symmetry line 
p — p x = Py', (b) the values of the electron Green's function << a prT \ a pa >> residues Z p ^ which correspond to 
the poles E^ (p) . Different line types represent four bands and corresponding residues. 

FIG. 4. The same as in Fig.3, but for truncated basis of three operators a p ,b p ,c p 
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